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Abstract 

We consider the problem of wavenumber selection for fully nonlinear, small-wavelength 
Gortler vortices in a curved channel flow. These type of Gortler vortices were first considered 
by Hall & Lakin (1988) for an external boundary layer flow. They proved particularly amenable 
to asymptotic description, it was possible to consider vortices large enough so that the mean 
flow correction driven by them is as large as the basic state, and this prompted us to consider 
them in a curved channel flow as an initial application of the phase-equation approach to 
Gortler vortices. This involves the assumption that the phase variable of these Gortler vortices 
varies on slow spanwise and time scales, then an analysis of both inside and outside the core 
region, to which vortex activity is restricted, leads to a system of partial differential equations 
which we can solve numerically for the wavenumber. We consider in particular the effect on 
the wavenumber of the outer channel wall varying on the same slow spanwise scale as the phase 
variable. 


^his research was partially supported by the National Aeronautics and Space Administration under NASA 
Contract No. NAS1-19480 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23681-0001. 
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1 Introduction 


Much innovative work was done during the 1980’s on wavenumber selection, particularly 
for convection problems. Attention turned towards selection mechanisms that relied on the 
global physics and geometry of the system. The method developed is now referred to as the 
‘phase-equation’ approach. Much of this work (see for instance Kramer et. al. 1982, Cross & 
Newell 1984, Buell & Catton 1986 and the review article by Newell, Passot & Lega 1993) was 
concerned with Rayleigh-Benard convection problems although the ideas behind it had been 
encountered earlier in the context of travelling wave instabilities (Whitham 1974 and Howard 
& Kopell 1977). 

Rayleigh-Benard convection is a classical example for pattern forming transitions in nonequi- 
librium systems. The spatially uniform conducting state becomes unstable to spatially periodic 
time-independent rolls. The stability of these rolls had previously been studied by Galerkin 
techniques which suggested that for a range of Rayleigh numbers (dependent on the Prandtl 
number and wavenumber) these straight rolls are stable. However, in experiments large enough 
to contain many rolls these solutions were usually not seen. Instead more complicated patterns 
are common with curved rolls, roll dislocations and superimposed rolls present and sometimes 
time independent states are not even reached. The presence of these complicated patterns was 
attributed to factors such as the existence of orientational degeneracy, the observation that 
rolls tend to align themselves normal to lateral boundaries and the fact that there is a band of 
stable wavenumbers. 

The ‘phase-equation’ approach did much to explain and predict the complicated patterns 
that are seen in experiments. For example, the dislocation of convection rolls is now fairly well 
understood using the phase equation approach. The phase-equation approach assumes that 
the phase is a function of slow (global) variables. All quantities are then expanded in terms of 
a suitable small parameter and substituted into the governing equations. At leading order the 
unmodulated equations of motion are recovered but at next order a linearised inhomogeneous 
form of the leading order problem is obtained. In order for this inhomogeneous problem to have 
a non— trivial solution a solvability condition must be satisfied and it is this solvability condition 
which gives an equation for the phase. 

Recently these ‘phase-equation’ ideas have begun to be applied to boundary layer prob- 
lems. Hall (1994) has used the phase-equation methods to consider large amplitude Tollmien- 
Schlichting waves in boundary layer flows. He considered asymptotically both large Reynolds 
number flows and finite Reynolds number flows (for the asymptotic suction profile). The re- 
sults suggest that for both large and finite Reynolds numbers a uniform wavetrain of Tollmien- 
Schlichting waves will break down with either a singularity or a shock developing after a finite 
time. Here we shall use these phase-equation ideas to consider wavenumber selection for small- 
wavelength Gortler vortices in a curved channel flow. 

Hall (1982) considered weakly nonlinear small-wavelength Gortler vortices, in an external 
boundary layer flow, and showed that this nonlinear interaction is not described by the Stuart- 
Watson approach. Hall & Lakin (1988) followed up this work by looking at the fully nonlinear 
interaction and they found that in the region of vortex activity the boundary layer flow is being 
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forced by the vortex which is in turn driven by the boundary layer. Surprisingly they also 
found that an asymptotic description could be obtained for this nonlinear interaction when the 
vortices are of sufficient size that the mean flow correction driven by them is as large as the 
basic state. This amenability to asymptotic description was one of the main reasons for making 
an initial phase-equation study on this type of Gortler vortices. 

In the next section we shall consider the formulation of our problem. iFrom the previous 
work of Hall k Lakin (1988) and Bassom (1989) we know that the vortices will be contained 
within a core region. In §3 we consider asymptotically the flow within this core region. In 
§4 we match the solutions from inside and outside the core region. In the process we derive 
a system of partial differential equations involving the wavenumber. Then in §5 we consider 
numerical solutions of the equations we have derived for a variety of different geometries for 
the outer channel wall. Finally in §6 we draw some conclusions from our numerical results and 
consider how far we have progressed in providing a theoretical explanation for the non-uniform 
patterns of vortices frequently observed in experiments (see Swearingen k Blackwelder 1987 
and references therein) and also in computational studies (Guo k Finlay 1994). 

2 Formulation of the Governing Equations 


In this work we consider the flow of an incompressible viscous fluid of density p and kinematic 
viscosity u in a curved channel. The walls of the channel, with respect to the usual cylindrical 
polar co-ordinates (r', 6' , z'), are given by r = R and r = R + dq(z) where d is a typical gap 
width (see the next section for a discussion of the slow spanwise scale, z = e*z, on which the 
function q varies). The non-dimensional co-ordinates ( x,y,z ) are defined by 


i m' 

X ~ Re d 



(2.1a, 6, c) 


where the Reynolds number Re = V a djv and V a is the maximum azimuthal flow velocity. 

The basic velocity and pressure fields for this small gap limit (6 = dj R <C 1) take the form 


(u,v,w) = V a (u,Re 1 v,Re 1 w) + 0(6V a ) , P = P 


( 2 . 2 ) 


and we confine our attention to the limit Re -* oo with the Gortler number G defined by 


, V 2 <P 


held fixed. We find that the Navier-Stokes equations for this problem are 


du dv dw _ 
dx dy dz 


d 2 u d 2 u du dp du du du 

dy 2 ^ dz 2 dt dx U dx + V dy W dz 


(2.3) 


(2.4a) 


(2.46) 
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(2.4c) 


d 2 v d 2 v dv 1 


dp 


_ _ , _ _ _ _ ‘n 2 _ _ 

O ^ yjli r\ 


dv dv 


dy 2 dz 2 dt 2 


u 


+ u— + u> 


dy dx dy 


dv 

dz 


d 2 w d 2 w dw dp dw dw dw 

+ -7TT - -«7 - a“ = U— + V— + W— 

dt dz dx dy dz 


dy' 2 ' dz 2 ~ a.. • ~ n_ • (2.4c?) 

We choose to concentrate on vortices with small, 0(e), wavelengths. Previous work has 
shown that such vortices exist at large, 0(e -4 ), Gortler numbers. So we expand the Gortler 
number for this flow as 

G = e -4 (G 0 + tG\ + • • •) , (2.5) 

and the mean flow quantities as 

(u, V, w ) = (u(y, 2, t), e^viy, z, t), e~?w(y, z, t)) 

P = -2 F(z, i)x + e~ 4 p(y, 5, t ) . (2.6) 

Again, see the next section for a discussion of the slow timescale, t = e 3 t . We find that vortex 
activity is restricted to a core region and outside of this core substitution of (2.6) into the 
governing equations, (2.4a-d), gives 


subject to the boundary conditions 


dv dw 

dy = ~~di ’ 

(2.7a) 

d 2 u 

dy 2 ~ ~ 

(2.74) 

<£|c£ 

II 

IS 

>— * I CN 
1 

(2.7c) 

II 

13 ^ 

CO <^) 

(2.7<f) 

w — 0 at y = 0, q . 

(2.8) 


3 Derivation of the Nonlinear Equations for Gortler 
Vortices in Curved Channel Flow 


We now develop an asymptotic solution of the governing equations valid in the presence 
of vortices which have small wavelengths. The previous work of Hall(1982), Hall k Lakin 
(1988) and Bassom (1989) suggests that the vortices will be contained within a core region. 
The flowfield is therfore split up as shown in Figure 1 with the vortex activity confined to the 
region between y\ and y 2 - The layers denoted by regions IIa,b are transition layers, of width 
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0(e f), required to smooth out the algebraically decaying vortices in region I. In region I the 
appropriate expansions of u, v, w and p are 


u = uo + e 2 u-[ + eu 2 H b 


V = e 3 Vo + ^Vi + €*V2-\ b 


w = e 2 tn 0 +u;i + e 2 u; 2 d b 


p = -2Fx+r 4 p 0 +e~^p,+e 3 p 2 + ■ ■ ■ 


e{U 0 E + c.c. ) + d(U l E + c.c.) + e 2 {U 2 E + c.c.) + 

e~ 1 (V 0 E+c.c.) + e~^(ViE+c.c.) + (V 2 E+c.c.)-\ +• 

(W 0 E-\-c.c.)+e^(WiE+ C . C .)+e(W 2 E+c.c.)-\ +• 


+ 


e~ l {P 0 E + c.c.) + C*{P X E + c.c.) + ( P 2 E + c.c.) + - - • 


+ 


, (3.1a) 


, (3.16) 
, (3.1c) 


(3.1 d) 


Mean flow terms (apart from the pressure gradient F) are denoted by barred small letters and 
vortex terms by capital letters. We assume that the mean flow quantities and disturbances are 
all independent of x and the notation c.c. denotes the complex conjugate, whilst the exponential 


E is given by 

E = expj 4-0(2,?))) 


(3.2) 


where 

z = th , t = e 3 t , (3.3) 

and 6 is a phase variable. The slow z-scale is determinded by our requirement that djdz ~ e -1 
and also by balancing terms in the mean flow z-momentum equation. In this paper we have, 
for convenience, fixed our attention on z = (az although we could have considered z = t k z 
where 3 < k < 4. Hall (1985) showed, in contrast to the two-dimensional case, that the Gortler 
disturbances applied to a three-dimensional flow are necessarily time dependent and hence 
the inclusion of this slow time-scale which is necessary to provide a term for the consistency 
condition at the second order of the fundamental equations (equation 3.13a). We write 0 = 
0{z,t) so that d/dz -> ike~ l d/dQ and d/dt -► -iCle~^djdQ where the wavenumber k = Q- z 
and the frequency ft = —9- t . The introduction of the wavenumber and frequency imply that we 
must have a phase conservation equation, namely 


dk dn n 

¥+a! = 0 


(3.4) 


Strictly, higher harmonics are also present in the expansions (3.1), but they do not affect our 
calculation at the orders we shall be concerned with. 

Substituting the expansions (3.1) into the governing equations we find that the leading order 
fundamental equations give 

Voy + ikWg = 0 , (3.5 a) 


k 2 U 0 + uoyVo = 0 


(3.56) 


k 2 V 0 + GqUqUq — 0 


(3.5c) 
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k 2 W 0 + ikP 0 = 0 , (3.5 d) 

and consistency of (3.5b) and (3.5c) implies that throughout the core 

G oUo^oy — k , (3.6a) 

from which we can see that 

+ a (5, 

ol 

where a is an, as yet undeterminded, arbitrary function of z and t. At the leading order of the 
mean flow we find that the governing equations give 


_ (2 k 4 y 

u 0 = 


*))* 


(3.66) 


Voy + Wqz — 0 , 

u 0 yy T 2 F ~ V 0 Uo y + V O^Oy + ikUoWo — i ku qWq , 
1 


(3.7a) 
(3.76) 

2 GoUo + p Oy = 0 , (3.7c) 

0 = V 0 W 0y + VoWoy , (3.7d) 

where we have used a bar on the capital (vortex) terms to denote the complex conjugate of this 
quantity. Using (3.5) to substitue in for Uq and Wo we find that (3.7b,d) give 

u 0yy + 2F = ~k^~dy ( u ° 1 '^’ 1 ^ 0 ) ’ 


rs / 

0=—\VoVoy-VoVo y 

and these two equations define the vortex function Vo- /,From (3.7c) we have 

Po y = -\(M 4 y + a) , 

and after integration with respect to y we have 

Po = -\{k A y 2 + ay + c(z,t)) . 

The size of the vortex function Vo is now determinded by integrating (3.8a) to give 

u 0y + 2 F y + B(z) = — -^uoy|Vo| 2 • 


(3.8a) 

(3.86) 

(3.9a) 

(3.96) 

(3.10) 


The function |Vo| 2 cannot be negative so yi and y 2 , which determine the edges of region I, must 
satisfy (3.10) with Vo = 0. If we then eliminate B(z) from the resulting two equations we obtain 
the condition 

k 4 k 4 

T +2Fi/ 1 = - T — ~+2Fy 2 . (3.11) 


Go{2k 4 y l + a) 2 


Go(2k 4 y 2 + a )2 
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This equation is obviously not sufficient to determine a, y\ and j /2 so we are not yet able to 
determine the location of the layers IIa,b. 

At second order the fundamental equations give 


V, y + *W, =0 , (3.12a) 

k 2 U 1 + u 0 yVi = -UiyVo - i(—0 + w 0 k)U 0 , (3.126) 

k 2 V\ -f- GqUqU\ = — GqU\Uq — t( — ft + U7ofc)Vo , (3.12c) 

k 2 W , + ikP \ = -i{-0 + w 0 k)W 0 , (3.12 d) 

and consistency of (3.12b) and (3.12c) requires firstly that 

— O + wok = 0 , (3.13a) 


and secondly that 


d_ 

dy 


(u 0 ui) — 0 


(3.136) 


Equation (3.13b), as we expect, defines the second order component of the streamwise mean 
flow whilst (3.13a) provides an equation for the leading order component of the spanwise mean 
flow. Note also that, since ft and k are not functions of y, equation (3.13a) implies that w 0y = 0. 

We now need to determine the positions of the transition layers, IIa,b. The thickness of 
these layers is determinded by a balance between diffusion across the layers and convection in 
the streamwise direction. This balance shows the layers to be O(es) in depth and so we define 
in region Ha the 0(1) variable 

^ = (y - 1 / 2 ) 


2 

63 


We are led to make the expansions 


u = (uqo + e 3 U 01 + c3Uq 2 + •••) + c 5 (^io -fi Gun + e 3 u 12 + ••■) + ••• 


(3.14) 


+ {e i >(EU l0 + c.c.) + e 2 {EUn + c.c.) + •••} + ••• , (3.15a) 

v = e 3 (noo 4 - e 3 Uoi + G V02 + **•) + e 2 (t>io + + c 3 Ui2 + •••) + --- 

+ {e->(EV 10 + c.c.) + (EV U + c.c.) + •■■} + ••■ , (3.156) 

W = 6 _ 2 (WOO + t 3 ^01 + t*W 02 + *•’) + (^10 + e*W U + £ 3 W\2 + •••) + -■■ 

+ {e"3 (EW l0 + c.c.) + r^(EWn+c.c.) + ■■■} + ■■■ , (3.15c) 

p = -2F:r+e -4 (p 00 +e3p 01 +e3p 02 -) )+e~ 2 (Pio+ e5 Pu+ e3 Pi23 )+' ' ' 

+ {e 3 [EP\o + c.c.) + e *(EPu + c.c.) + ••■} + ••• , (3.15d) 

in this layer because we can see from (3.10) that |Vo| 2 ~ P 2 — V when y — * 3 / 2 — - 
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Considering the solution of the governing equations in this layer we find that the fundamental 
x- and y-momentum equations give at 0(e - 3) and 0(e _ 3) respectively 

k 2 U w + Vio q = 0 , (3.16a) 

k 2 V\o + GqUqqU^o = 0 , (3.166) 

whilst at 0(1) and 0(e -2 ) of the same equations we obtain 

1,2j T , T/ ^ 01 _ PUlO T/ du 0 2 /o 1-7 \ 

k U u + - V 10 — , (3.17a) 

— d 2 V 

fc 2 Vii + GqUooUw = 2 2 Goaoif/io • (3.176) 

Equations (3.16a,b) are always consistent but for (3.17a,b) to be consistent we require that 


, ,_4 U 01 t/ , n = #U02 

3£ — h k = — Vio + OoWoo _ "x t Ko = 0 


(3.18) 


Looking at the mean flow x-momentum equation we find that at 0(e a), 0(e 3 ) and 0(1) 


d 2 u 00 _ n 

ae ’ 


^ 2 U 0 l 


T = ° . 


9 2 «02 , or_i/ WlO , T7 dU 10 , • j jj 7 rr 77 

\- 2 F — Vio— — h V 10 — 1“ ikW iqU\o — ikWwU 10 


The solutions of (3.19a) and (3.19b) which match with the solutions in the 


core are 


= ( 2k 4 y 2 + a)i 

Uoo = I 


(3.19a) 


(3.196) 

(3.19c) 


(3.20a) 


Uoi — 1 ) 

Go{2k 4 y 2 + a)3 


(3.206) 


whilst we find from (3.19c), after using the fundamental continuity and x-momentum equations 
to substitute in for U 10 and W/i 0 , that 


9^02 = _ 2 P d|Vjol 2 

Gi(2k*y 2 + a)i 


(3.21a) 
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and integrating this with respect to £ gives 


= -m 


2k 2 \V v 


Go(2k 4 y 2 + a 


— + f(z) , 


(3.216) 


where f is an arbitrary function of z and t which we do not need to know explicitly for our 
purposes. If we now substitute (3.21b) back into equation (3.18) we find that 


-glgr ■*" ts(iMo = + g{z,i)v. k 


where 


S{z,t) = 


k r 6 _ 2FG|(2fc 4 y 2 + a)i 

3(2k 4 y 2 + a) 3 k 2 


- Gl{2k A y 2 + a)h 

9W) = 3^ J 


(3.22) 


(3.23a) 


(3.236) 


Equation (3.22) is a particular form of the second Painleve transcendent and has been shown 
by Hastings & McLeod (1980) to have a solution such that 


as £ — > —00 


0 as £ 


(3.24) 


It follows that in layer I la the fundamental terms decay to zero so that the finite amplitude 
Gortler vortex is trapped above region Ilia. We note that a similar analysis for the higher 
harmonics shows that these functions also decay exponentially to zero in Ha. However, the 
mean flow is virtually unaltered by the presence of Ila, thus the early terms in the expansions 
of the mean flow quantities in region Ila are simply obtained by expanding the mean flow in 
region I in terms of £. This means that the mean flow in region Ilia must to zeroth order when 
y —> t/ 2 + have u, u y , v, w, w y and p defined by the coreflow solution evaluated with y - y 2 . 

An identical analysis can be carried out in region lib with similar results. Hence in both 
region Ilia and Illb there is only a mean velocity field. Therefore in the regions (0,j/i) and 
(y 2 ,q) the equations (2.7a-d) hold and we must solve them subject to the boundary conditions 
(2.8) and the matching conditions 


u{yi) = 


(2k 4 yi + a) 2 


Uj ,(yi) = 


Gq (2 k 4 yi + a): 


w(yi) -Q/k , WyiVi) = 0 , 

v(yi) = —yi(Q/k)z + 6 , 


p(yi) = + (iy, + c) 


t = 1,2 


(3.25a, 6) 

(3.25c, d) 
(3.25e) 

(3.25/) 


along with the relationship between the positions y\ and y 2 (equation 3.11). 
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4 Derivation of the Vortex Wavenumber Evolution Equa- 
tion 


Above the core, in region Illb, equation (2.7b) gives 

d 2 u 


dy 2 


= -2 F(z,t) , 


(4.1) 


and integrating this twice with respect to y gives, after we have applied the boundary condition 
u = 0 at y = 0, 

u = y(A - yF) , (4.2) 

where A is an as yet undeterminded arbitrary function of z and t. Matching u and u y at y x 
(see equations 3.25a,b) implies that 


k 4 = Goyi(A — y x F)(A — 2y x F) , 
a = -G 0 y 2 1 (A-y 1 F)(A-3y l F) . 

Below the core equation (2.7b) again holds and integrating twice gives 

u = - Fy 2 + ay + b , 
whilst matching u and u y at y 2 we find 

yl{A - y x F) 2 (A - 2 y x F) 


a — 


b = -Fy] - 


[2y 2 (A-2y 1 F)-y 1 (A-3y 1 F))2 

y 2 y\{A-y l F) 1 2{A-2y l F) 


t + 2Fy 2 


[2y 2 (A-2y 1 F)-y l {A-Zy 1 F)]$ 

+yt(A - yiF)s[2y 2 {A -2y 1 F)-y 1 (A- SyiF)]! . 

We must satisfy the boundary condition u = 0 at y = q and this implies that 

0 = — F y 2 -(- ay b , 

and upon substituting in for a and b we have 

o = - F (q - „ 2 ) 2 + 

{2y 2 {A-2 yi F)- m (A-3 yi F)]i 

+y'l(A - y,F)i[2y 2 (A - 2 y,F) - y,(A - %,F)]5 


(43) 

(4.4) 


(4.5) 


(4.6) 

(4.7) 


(4.8) 
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^From our working inside the core we have equation (3.11), the relationship between the po- 
sitions 2/1 and y 2 , and after substituting in for k 4 and a, using equations (4.3) and (4.4), this 
becomes , 

(A - 2 y,F) = i lHA-ViF)HA-1viF) ( 4 . 9 ) 

{2y 2 (A-2y,F)-yi(A-3y,F)\l 

We now have three equations (4.3, 4.8 and 4.9) with which we can determine A , j/i and y 2 in 
terms of k and F (assuming we know G 0 and q ). 

Above the core the y-momentum equation is 


dp l r _ 2 

Ty 2 G ° U ’ 


and after substituting for u and integrating we find that 

p = -IGoYtF'v* - l -AFy l + X -A 2 y 3 \ +Q('z,t) , 


2 "5 


(4.10) 


(4.11) 


where Q is an arbitrary function of z and t. Matching with the pressure in the core at t/i we 
find that the constant c is given by 


c = Go - ^AFyf + -2Q 


(4.12) 


Below the core we find that 


p= -i^o \F‘ 2 y 5 -\aFy 4 + \{a 2 -2bF)y z + aby 2 + b 2 y + c(z,t) , (4.13) 

Z [o Z 6 

and matching with the pressure in the core at j /2 we get 

c = iGofiF 2 v|-iaFv‘+i(a 2 -25F)l(!+«5!li+6 ! <(2-5F 2 ^+|2lF! ) ;-i2l 2 yf 


2 5 


-y\yi{A - y\F){y 2 {A - 2y x F) -yi(A-3yiF)} +Q . (4.14) 

The z-momentum equation in region Illb is (2.7d) and integrating twice with respect to y 
and using the boundary condition w = 0 at y — 0 we obtain 

W = ^{F 2 };y 7 - -{AF}- Z y & + —{A 2 }^ 5 + -Qzy 2 + C(z,t)y , (4.15) 

where C is another arbitrary function of z and t. We know that the position of the upper edge 
of the core occurs where 8w/dy — 0 and this produces the condition 

0 =-5Go[4 { F 2 }«f-4 { 4F}^ + i { X 2} , y ;]+(? i y 1 +C . (4.16) 
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From matching w at this upper edge of the core we obtain 


1 

I “ ~2° 


2l0 + 


+ yQzy* + C3/1 • 


Integrating (2.7d) in region Ilia twice gives 


= _ 5 Go 


!} » 7 - + - L {« 2 - 2 bFW + + i { 5 J },-!, 3 


+\c- z y 2 + d(z,t)y + e(z,t) , 


and matching dw/dy at 7/2 gives 


3 =- G ° 


^{Fn^!-^{aC}«l+^{3 2 -26F}^+i{a6}«I + i{S 2 }rt ! 2 


To 


12 


1 

3 


whilst matching w gives 


- « 

'“I + 2 Go 


- ^{Sf}i»| + ^{<f 2 - ibF}iVl 


(4.17) 


(4.18) 


cjt /2 , ( 4 . 19 ) 




!_ 2 -* 

- 2 C * y * ~ dy2 


(4.20) 


We must satisfy the boundary condition w = 0 at y = q and from equation (4.18) this implies 
that 


0 = —Go 




+r^{ab};q 4 + b 2 } z q 3 


1_ 2 -= 

+ -c z q + dq + e . 


(4.21) 


Finally we need to consider the continuity equation. Above the core we find, after integrating 
(2.7a) and applying the boundary condition v = 0 at y = 0, that 


H g ° 




- IfeV - \c-tf , ( 4 . 22 ) 


and matching v at y i implies that b is determinded by 


6=i Go 


1680 ^ F ^ hV ' 420 + 360 


--^QzzVi — 2Czy\+{£ljk}z y i . (4.23) 


In region Ilia we again have (2.7a) and after substituting in and integrating once we find that 


”=2 G » 


1680 {F 2 }zzV B 420^ aF ^' i2/7+ 360^ 2 2bF }™ y6 + 60 ^°^ 55y5+ 24 ^ 
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(4.24) 


-jC.'iJ 3 - 7, d iy 2 - Ci'J + f(iJ) 
whilst matching with the vertical velocity in the core at y 2 determines /; 


7= \ a ° 




-2 G ° 


1680 }~ z ' zV2 420^ aF ^' zzy2 


+ ^{a 2 - 2 bF}- z - z yt + ^{ab};zV 5 2 + -^{b 2 }t-zyl 


+ {£l/k}z(yi j/2) qQzzV\ 


1_ . 1-= 2 

~^zzV2 n d zU2 “I” 2 

0 Z 


(4.25) 


The final boundary condition we have to satisfy from (2.8) is v — 0 at y = q and this implies 
that 


0 = 3 Go 


16^ <f2}a “430 < “ F, " ,7+ 3k {aJ_25F) " ,6+ S {ai ' } “ ,S+ S !i,2}lM4 


-4ch? 3 - l-dsq 2 -e t q + f . (4.26) 

0 Z 

We now have seven equations (4.3, 4.8-9, 4.16-17, 4.21 and 4.26) which relate the eight unknowns 
(A, F, C, Qz,Vi, 2 / 2 , and k). If we can determine a solution at an initial time location we can 
then use the phase conservation equation 


dk dfl 

dl + as =0 


(4.27) 


to march the solution forwards in time. 


5 Numerical Solutions 


In this section we explain the numerical method used to obtain solutions of the system of 
equations derived in the last section and we consider results for a variety of different outer 
channel walls and initial wavenumber distributions, k(z,t = 0). 

If we rescale A, y\ and y 2 by writing 

A = qFA , j/i = qyi , y 2 = qh , (5.1) 

we find that equations (4. 3, 8, 9) become 

2S,) . (5-2) 

G 0 q 3 F 2 

Q = -d - y 2 f + (1 ~ Mfh* - ~ 

[2 y 2 (A - 2t/i) - yi{A - 3yi)> 
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(5.3) 


+Vi (A - yi ) 2 [2 y 2 (A - 2y, ) - y x (A - Zy x )] a , 

yM-yi)M-2yi) 

[2y 2 (A - 2ft) - ft(A - 3t/! )] 2 


(5.4) 


We can introduce the variable 77 = k 4 / ( Goq 3 F 2 ) and then solve for A, ft, ft and their derivatives 
numerically. We do this by taking an initial guess at ft (for a given 77 ) so that (5.2) is a quadratic 
in A. After solving this we then solve the square of (5.4) which is a cubic for ft and substitute 
our answers into (5.3) to check if we need to iterate on our choice of ft . Obviously there are six 
possible sets of answers but we find that only one set is permissible (both ft and ft real with 
0 < ft, ft < 1 and ft > ft). We plot A, ft and ft along with their first and second derivatives 
(as these will be needed for our further calculations) in Figures 2-7. Numerically we find that 
a solution exists for 77 greater than zero and less than approximately 0.104 . The lower limit 
corresponds to the core occupying the entire channel and the upper limit occurs when ft and 
ft coalesce and the core region disappears. However, we note that there may be some error in 
the value of this upper limit as the solution becomes extremely difficult to follow numerically 
in this region because of the rapid changes in the derivatives. 

If we know the distribution of k at an initial time and also F and its first derivative at 
some z location we can use equation (4.21) to determine the unknown Q 5 . Upon substitution 
equation (4.21) becomes 


Qs = 


[(9 - y-iY - y\\ [ 2 L210 


G 0 


1 ;<r 2 } i (s<f[7?-6!/2]-9 7 )-4j{a7 ; ’}i(y5[69-5!/j)-}' i ) 


+^{(5 2 - MFMvUSi - 4y 2 ] - 9 s ) + - 3^] - 9 4 ) + - 2y»] - ? 3 ) 


l [q 2 12^ AF ^ !V ' + 


(5.5) 


where for convenience we have introduced 


c — c — Q . (5.6) 

We can then determine the unknown C from (4.16) and the frequency ff from (4.17). Finally 
after lengthy substitution we can determine the second derivative of F from the final boundary 
condition, (4.26), and thus obtain a full solution of our equations at an initial time location. 
This solution is then marched forward in time by using the phase conservation equation, (4.27). 
We note that (5.5) implies that a singular solution will occur for Qz when y\ — (q — y 2 ) but, 
this only happens when 77 = 0 (see Figure 2 ) and so will not concern us. 

We now consider some specific functions for the outer wall, q , and determine in particular 
the effect they have on the vortex wavenumber. For all of these calculations we shall take 
Go = 1 and assume that F(z = 0) = 1 , F^{z = 0) = 0. Obviously there is a virtually infinite 
choice of outer walls and initial wavenumber distributions and our numerical work here only 
covers a few of these. The particular cases we highlight here are relatively simple ones and 
have been chosen in an effort to display the effects that individual features have. We have 
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concentrated on variations in the outer wall as opposed to variations in the distribution of 
the wavenumber at t = 0 (which we have taken to be constant in all of these cases) because 
the wall geometries we consider cause the wavenumber distribution to become non-constant 
anyway. We were also particularly interested in the effect that the geometry of the channel has 
on the observed wavenumber. 

The first case we considered was that where there is a slight increase in the width of the 
channel. The function 9 was defined as 

9 =1.15 , z <10.0 and q = 1.2 - 0.05e- (z '- lo)3 , 5 > 10.0 , (5.7) 

and is shown in Figure 8 along with the extent of the core region for an initial distribution of 
the wavenumber Jc(z, t = 0) = 0.5 . In Figure 9 we plot the value of the wavenumber at time 
intervals of t = 125, 000. The expansion of the channel causes the wavenumber to increase before 
the expansion and decrease after it. As time passes the effect of the change in the channel width 
on the wavenumber is felt both further upstream and downstream in the spanwise direction but 
we note that the wavenumber changes very slowly in time. The changes in the wavenumber 
are insufficient to produce a change in the core position that would be graphically noticeable 
on Figure 8. These results plus further unplotted ones and a consideration of dCl/dz values 
suggest that the wavenumber eventually settles down to a different steady state. 

The second case we considered was that of a constriction of the channel with the function 
q now given by 

9 = 1.1 , z < 10.0 and q = 1.0 + 0.1e-< 2 ~- 10 ) 3 , z > 10.0 . (5.8) 

See Figure 10 for a plot of this and the position of the core, where we have again choosen an 
initial distribution of the wavenumber k(z,t = 0) = 0.5. In Figure 11 we show the value of the 
wavenumber at t = 125,000 and t = 250,000 for this case. A constriction of the channel has 
the opposite effect to an expansion with the wavenumber decreasing before the change in the 
channel width and increasing afterwards. The other point to be noted from this case is that 
a larger change in the channel width (0.1 as opposed to 0.05 for the expansion of the channel 
that we looked at) has produced a larger and quicker change in the wavenumber. 

A more interesting case to consider is that where the outer wall of the channel is a periodic 
function. For this we have considered 

9=1.8 , z < 10.0 and z > 10.0 + 4tt , 

9 = 1.775 + 0.05 cos(z — 10) — 0.025 cos 2 (i — 10) , 10.0 < z < 10.0 + 47r , (5.9) 

so that the first and second derivatives of q are continuous throughout the range of 5 that we 
consider. Figure 12 shows the channel and core position for this wall function and, again, a 
distribution of k = 0.5 at t = 0. We note that the top edge of the core is sufficiently close to the 
wall to not be visible in this graph and also that the position of the core continues to change 
for values of z considerably greater than that at which changes to the channel wall stop. The 
value of j/ 2 > the lower edge of the core region, eventually settles at about 0.931. In Figure 13 we 
have plotted the wavenumber at i = 4, 000 and 7 = 8, 000. The wavenumber is again changing 
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slowly in time and we also note that in the spanwise direction it soon returns to its initial value 
after the channel wall stops altering. These results bring up the interesting possibility of some 
complicated wall variation having a small effect on the vortex wavenumber itself (and only close 
to the position of the variations) but a much larger influence on the vortices through its effect 
on the position of the core region, maybe even leading to the core region being pinched out of 
existence at a location some distance downstream (in the spanwise direction) from the changes 
in the channel geometry. 

Further investigations show that the disappearance of vortex activity altogether through 
the narrowing of the core region is possible but this will probably occur close to the area where 
the wall alters. If we change the initial distribution of the wavenumber to k = 0.70 we find (see 
Figure 14) that although we end up with a narrower core region the changes in the position of 
this core region stop not long after the area of wall variation is passed. Increasing the initial 
value of the wavenumber further to k = 0.7092 results in the core region of vortex activity being 
extinguished, Figure 15, but this occurs just after the region where the outer wall is periodic. 


6 Conclusions 


In this paper we have considered the problem of wavenumber selection for small-wavelength 
Gortler vortices in a curved channel flow. Following an analysis similar to that of Hall &; Lakin 
(1988) but with modifications for the flow being in a curved channel and the vortex wavenumber 
being dependent on slow spanwise and time scales we have obtained a system of equations which 
can be solved numerically for the vortex wavenumber. Our results, contained in the previous 
section, concentrate on the problem where the outer wall of the curved channel varies on the 
same slow spanwise scale as the phase variable of the vortices. The reasoning behind this was a 
desire to determine what effect wall geometries could have on the patterns of vortices observed. 
Alternatively we could have considered the evolution of a system of vortices that started, for 
our problem, with a non-uniform spanwise wavenumber. 

The work in this paper has shown that it is possible to use phase-equation methods to 
track changes in the wavenumber of Gortler vortices. It would be desirable to extend this to 
the problem of an external flow where the vortex wavenumber would be dependent on slow 
time, streamwise and spanwise scales. In particular this would complement other recent studies 
on Gortler vortices in three-dimensional boundary layers which have been prompted by the 
development of laminar flow control airfoils with significant areas of concave curvature on the 
underside of the wing near to the leading edge. Unfortunately it is not possible to consider the 
current problem in an external boundary layer because the structure of equations (3.1a-d) will 
not suffice when the phase variable of the vortices is a function of the streamwise variable. So 
in order to make such an extension it will be necessary to find an alternative system of Gortler 
vortex equations that is ameniable to asymptotic investigation. We note however, that this may 
be difficult and even for the problem considered in this paper, which is particularly receptive 
to an asymptotic analysis, we end up with a complicated system of equations; the solution of 
which requires quite lengthy numerical computations. 
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The numerical results we have obtained suggest that the vortex wavenumber will respond 
to variations in the geometry of the outer wall of the channel with small changes occuring 
over quite large lengths of time. These small changes in the wavenumber cause only minor 
changes in the other flow quantities as time progresses. However, more dramatic changes in 
the spanwise direction for the flow outside the core (caused by the variations in the wall) can 
it appears influence the vortices. It seems that for certain initial conditions and wall variations 
(see Figure 15) the core region of vortex activity can even be pinched out of existence. 

Finally we compare our results with those of experimental studies on Gortler vortex wave- 
length selection. The large size of our Gortler number makes direct comparison with experi- 
mental work difficult but it appears that this choice of Gortler number is required in order to 
make analytical progress. Our results (see Figures 9, II and 13) suggest that changes in the 
‘global’ geometry of the channel will produce different vortex patterns but will have very little 
effect on the average wavenumber (and hence wavelength). This is supported by the work of 
Swearingen & Blackwelder (1986); one of whose findings was that the average wavelength is 
independent of the test section spanwise width. 

The comparitively small changes to the vortex wavenumber which we have found in this 
paper show the need for a theoretical investigation of wavenumber selection at the onset of 
vortex instability in order to explain the large differences in the vortex wavelength which have 
been obtained in experimental work. It was observed by Tani (1962) and Tani & Sakagami 
(1964) that the vortex wavelength changed when a different experimental facility was used and 
it has been possible to alter the observed wavelength by the artificial means of the introduction 
of disturbances upstream of the onset of vortex activity (see Myose & Blackwelder 1991 for 
a discussion of the different experiments that have been conducted). However, as pointed 
out by Myose & Blackwelder (1991) the exact process of wavelength selection is not yet well- 
understood. Our results here confirm that the initial selection of the wavenumber, at the start 
of vortex activity, needs theoretical consideration. Unfortunately ‘phase-equation’ methods 
would appear to be of no use in describing the effects of disturbances that are periodic in 
the spanwise direction with spacing of the order of the wavelength of the Gortler vortices, as 
was the case for the experiments of Swearingen & Blackwelder (1986). Once this problem of 
wavenumber selection near the onset of vortex activity has been tackled the phase-equation 
method detailed in this paper will allow further changes in the wavenumber to be predicted. 
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Figure 1. A schematic picture of the different regions of the flow in the curved channel. 
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Figure 4. The second derivatives of y\ and y 2 
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Figure 5. The function A. 



Figure 6. The first derivative of A. 
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Figure 7. The second derivative of A 
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Figure 8. The function q from equation (5.7). 
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Figure 9. The wavenumber, k, for the channel given by q from (5.7) 
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